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ABSTRACT 

The  ratio  cf  geostrophic  drag  coefficients  C/Ca  under 
various  stability  conditions  is  obtained  as  a  function  of 
Richardson  number.   A  dra(  coefficient  Cf  is  defined  in 
relation  to  the  mechanical  mixing  length,  and  its  ratio  to 
the  actual  drag  coefficient  C  is  related  to  the  Richardson 
number.   By  using  observational  recent  data  on  the  normal- 
ized logarithmic  wind  shear  of  Monin-Obukhov  one  can  obtain 
Cf/Ca. 

The  writer  is  deeply  indebted  tc  Dr.  P.  L.  Martin 
(Professor  of  Meteorology)  for  his  suggestion  of  the  topic 
and  his  continued  help  throughout  the  investigation  and 
during  the  preparation  of  this  paper. 
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1.   Introduction. 

One  of  the  important  problems  in  meteorolo,' y  is  to 
establish  relationships  Vetween  the  turbulent  characteris- 
tics of  the  lower  atmosphere  and  the  larre-scale  synoptic 
parameters.   A  si  nif leant  turbulent  characteristic  is  the 
round  drag,  which  enables  one  to  specify  the  wind  profile, 
the  eddy  diffusivity  and  the  local  energy  dissipation  in 
the  surface  layer,  by  knowing  the  topof raphical  characteris- 
tics i.e.,  the  rovhness  parameter  of  a  particular  place. 

Lettau  [7J  had  made  an  attempt  to  relate  the  draf 
coefficients  with  the  static  stability  represented  by  a 
dimsnsicnless  parameter  called  Richardson  number.   His 
invest!  ations  were  based  on  1955  O'Neill,  Nebraska,  data. 
It  is  the  main  purpose  of  this  paper  to  verify  and  extend 
his  investigations  using-  1956  O'Neill  data.   It  is  important 
to  note  that  this  wo  rk    deals  with  a  relatively  smooth 
surface,  ti  erefore  the  "zero-plane  displacement"  which  is 
sometimes  used  in  the  adiabatic-wind  profile  nas  been 
neglected,  a  procedure  which  greatly  simplifies  the  com- 
putations. 


2.   The  Nature  of  the  Data. 

The  data  used  in  this  study  were  obtained  by  the 
Project  Frairie  Grass  flj  conducted  near  O'Neill,  Nebraska 
during  the  summer  1956.   Since  this  study  rakes  use  of  some 
of  the  published  data,  i.e.,  wind  speed,  temperature  and 
surface  maps,  it  is  necessary  to  describe  some  aspects  of 
the  instrumentation.   The  instrumentation  employed  was  that 
of  the  mobile  micrometeorological  station  of  the  Texas  A&M 
f roup.   The  station  has  a  slender  aluminum  mast  supporting 
six  anemometers  at  heights  of  8,  4,  2,  1,  0.5  and  0.25  meters 
A  similar  mast  supported  seven  temperature-measuring,  radia- 
tion-shielded copper-constantan  thermocouple  junctions  at 
heights  6,  4,  2,  1,  0.5,  0C25  and  0.125  meters.   In  addition 
the  published  report  of  Project  Prairie  Grass  includes 
sectional  sea  level  pressure  maps  at  times  nearlv  synoptic 
with  the  micrometeorolop ical  data.   These  maps  revealed  some 
small-scale  features  of  the  circulation  near  O'Neill  which 
did  not  appear  on  the  larger-scale  facsimile  maps0   The  micro- 
meteorological  data  oeriods  were  selected  to  be  nearly  simul- 
taneous with  the  available  maps.   At  tines,  linear  interpola- 
tion between  the  available  charts  was  used  for  computations 
of  geostrophic  wind  speeds. 

Regarding  the  errors  of  measurement,  two' types  of  tem- 
perature error  occur:  i.e.,  calibration  and  radiation  error0 
An  estimate  for  calibration  error  was  0.05OC.   The  estimated 
radiation  error  was  0ol°C  with  the  understanding  that,  in 
daytime  with  clear  skies,  and  low  windspeed,  all  measured 
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air  temperatures  overestimate  the  true  temperature.   At  night 
with  clear  skies  the  measured  temperatures  would  re  lower 
than  the  ambient  temperature.   For  windspeed,  measurements 
were  dependent  upon  different  technicues  of  calibrations 
according  to  the  windforce.   All  measurements  used  in  this 
paper  were  twenty-minute  means.   The  reader  is  referred  to 
[l]  for  more  details  concerning  the  instrumentation. 

5.   Background  Theory. 

Extensive  use  will  be  made  of  Prandtl's  turbulence 
model,  because  of  its  relative  simplicity,.   This  model  is 
described  in  many  standard  texts.   According  to  this  theory, 
the  followin  major  results  emerge: 


a  =ri& 


dz  <!> 


oVdu  \* 


where  JUL*  ,  OJ"* ,  are  the  tiirbulent  velocity  fluctuations  in 
the  direction  of  the  mean  wind  and  in  the  vertical,  respec- 
tively. ]L    is  the  eddy  mixin^  length,  while  a/  ia  the  root- 
mean-snuare  mixin?  len  th.   The  eddy  stress  X'  is  then 
considered  to  he  very  nearly  a  constant  X0  within  a  layer 
called  the  surface  layer. 

In  the  case  of  the  fully- rou^h,  neutral  surface  ]ayer, 
rtos^by  has  shown  that  As  Kl,  so  that  (2)  may  be  integrated 
to   ive 


LL*  =    -  U<  4 


where   UL^^vTc/p   Is  the  friction  velocity  in  a  neutral 
surface  layer,  Z0  is  the  ro^  bness  parameter  and  k=0.38 
is  the  Von  Karman  constant. 

In  the  non-neutral  surface  layer,  numerous  theoretical 
profiles  have  been  derived.   However  Deacon  has  proposed  the 
semi-empirical  relationship 

-(* 


(4) 


w 


do.     __  u* 

here   tine    friction   velocity  u.„,  sNi^/f      is    constant    in   the 


(*)' 


non-neutral  surface  layer.   The  parameter  6  is  a  decreasing 
function  of  Richardson  number.   Actually  Davidson  and 
T  arad  [4]  have  shown  that  \ G  - 1\  increases  somewhat  with 
hei  ht.   However  (4)  still  affords  a  useful  tool  as  an 
overall  or  hulk-relationship.   Two  useful  res;  Its  may  he 
derived  from  (4).   These  are: 

(i)  the  integrated  wind  profile: 

•-ft 

kd-fl) 

(ii)  the  value  JUL^  by  the  limiting  process  indicated 
in  eouatlon  (7)  below 


kd-fl)     \E*' 


(5) 


so   that 


1    f.       Ju76(M 


U*  —     K    w*^    -773 — \ 


(5) 


(7) 


Actually   the    lowest    level   at  which  wind    data  were   available 
was    Z =  25   cm    and    Z0    is    the    rouphness   ja  rameter  determined 
under   the   time-nearest   neutral    conditions,,   IX.^  may   therefore 
be   approximated   by    tl  e    finite-difference   expression 
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u  _    o.38aay 


Pol  low  in/,  a  procedure  similar  to  that  employed  by  David- 
son and  Farad  [4]  in  neutral  cases,  it  may  be  shown  that 


^i,- 


where 


A(W)  =  U.   800  +  ^  <_oo  +  ^  ^        (1Q) 


and 


all  =  (u«-u.1)>^-ua)+K-S)      (ID 

for  mean-wind  measurements  at  25,  50,  ICO,  200,  400,  800  cms„ 

4.   Geostrophic  Drag  Coef ficientse 

Lettau  l"7]  '  as  defined  the  '  eostrophic  drag  coefficient 
as  * 

C    =  a.  /v,  (12) 

For  neutral  stability,  the  corresponding  drag  coefficient 
is  defined  as 

Neutral  wind  profiles  at  O'Neill  may  be  obtained  by  finding 
those  cases  for  which  the  potential  temperature  profiles 
we  characterized  byd©/dzs0  in  the  layer  25  cm  -  200  cm0 
Altogether  15  such  profiles  were  found.   Assuming  a  logarith- 
mic wind  profile  can  be  fitted  to  the  wind  data,  U*^  of 
equation  (Z)   may  be  obtained  as 

UL ^U    -    V^  (14) 


By   analogy  with    flow    in   circular   conduits,    Lettau    [7] 
has    suggested   a   drag-coefficient   Ca   for  neutral  esses   of 
the   form 


Q.IOZj 


lo3(C«.Ro)-2>*Z, 


(15) 


where    do  may  be   called   the   "surface    Rossby   number"   defined 
as    the   non-dimensional    expression 


Vi=  4^ 


(16) 


*o   =V3/fZo      ,   V3.  ^   ^ 

The  constants  of  ecuation  (15)  were  obtained  by  Lettau  by 

least-square  methods  applied  to  a  wide  variety  of  neutral- 
wind  profiles  for  different  sites. 

Lettau  \_7j    goes  on  to  five  values  of  C/Ca  for  different 
values  iii]_oo»  the  Hichardson  number  at  100  cm,  for  the  1952 
O'Neill  data.  He  obtains  a  curve  as  shown  in  fig.  1„ 

1,3 


.COS    -'O-Os     -°°i     -**i    -°  Ol        o        O-Oi        ©.02       O-Cl       004    o-os~ 

■p" 

Fig,    1. 
Ratio   of  ceostrophic   drag   coefficient   for  non-adiabatic   and 
adiabatic    conditions    as   a    function   of    rtichardson   number  at 
Z=100    cm,    at   O'Neill,    Nebraska    1953    (After  Lettau  [7]  ) . 

It  must  be  pointed  out  that  Ca  here  represents  the  geostro- 
phic  drag  coefficient  for  any  surface  layer  of  drag  coeffi- 
cient   C,    if   this    layer  were    suddenly   to   be    converted    into 

5 


a  neutral  layer  without  changing  Vf/fz0. 

A  primary  objective  of  this  paper  was  to  test  eouation 
(15)  as  a  means  of  determining  Ca  based  on  the  1956  O'Neill 
data.   Also  application  of  non-neutral  orofile  data  to  eoua- 
tions  (7)  and  (12)  afforded  values  of  C  under  identical 
roufhness  conditions.   The  comparison  of  Lettau's  1953  graph 
of  C/Ca  versus  tfiioo  with  that  obtained  by  the  writer  at 
iii]_40  (to  be  more  exact,  the  Richardson  number  at  141 «,  4cm) 
based  on   the  1956  O'Neill  data  is  shown  in  fig.  3. 

If  one  knows  a  value  of  do   •=.   Vf/fz0,  it  is  possible  to 
solve  eouation  (15)  for  Ca  by  numerical  methods.   It  is 
shown  in  the  Appendix  that  a  numerical  solution  Ca  of  eoua- 
tion (15)  may  be  obtained  iteratively  from  the  eouation 

y-  (a  3026  ^  >;  +  V:  +  Cs^iy  -  z.ionflcj  "Ro  ) 


y>.  *  *  -  - 


»  +  y,- 


(17) 


where  vj_  =  i-th   iterant    to   0.2395/Ca.      The    first    estimate 
for  Ca  was   obtained    from   Lettau's    curve    [7]  ,    fig.    2, 


c.oso 


C..C-1&    . 


0..0'.0 


£?     So      5!?      to      ZT      7.0      j.s      8.0      8-s       $  0      f.»- 

Fif.    2S  ^Ro 

Geostrophic   drag    coefficient   under  neutral    condition   as    a 
function   of   surface    Rossby   number    (After  Lettau   [7]  ) 0 


The  chief  inaccuracy  in  this  approach  appears  to  be  uncer- 
tainty In  the  numerical  constants  appearing  in  eouation(15) , 
especially  at  low  geostrophic  wind  speeds.   Another  source 
of  error  is  the  necessity  in  some  cases,  of  interpolating 
Vg  tetween  two  maps  5  or  12  hours  apart  in  order  to  get  the 
estimated  value  of  Vg  at  the  time  of  microcieteorological 
data- run. 

A  second  approach  to  determine  the  relationship  between 
C  and  Ca  was  used.   Based  on  15  neutral  cases  at  O'Neill, 
Nebraska,  values  of  Ca  were  computed  using  equations  (13) 
and  (14).   Values  of  Ca  were  then  plotted  versus  log  Ro, 
and  the  resulting  best-fitting  curve  compared  to  Lettau's, 
which  is  replotted  on  the  same  graph.   The  comparison  will 
be  discussed  in  later  section. 

The  Richardson  number  has  been  used  very  extensively 
as  a  measure  of  stability  in  the  micro-meteorolory0   Its 
definition  is 

"     e(d&./*z 

This  number  is  proportional  to  the  ratio  of  work  done  against 
static  stability  to  the  work  done  by  eddy  stress. 

For  instruments  arranged  at  eoual  logarithmic  spacim, 
as  in  the  1955  and  1956  O'Neill  data,  a  convenient  finite- 
difference  form  of  Ri  at  Z=Z]_#2  nas  been  proposed  by 
Lettau  [s] 

fcfej,     *'-i?'.-*).^t  (19) 

This   value   is   understood    to    be    api  licable   at    the   f eometric 

mean  of   Z^   and    Z2,    namely   Zi.2  r  M  Z1Z2'.     @2   and     ^-o   in   the 
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Q(^/dzf  U8) 


1956  data  are  the  20-rninute  mean  potential  temperature  and 
windspeed  at  level  2,  and  an  analogous  meaning  applies  for 
symbols  with  subscript  one. 

Finally,  Monin  and  Obukhov  [VJhave  introduced  a  use- 
ful non-dimensional  parameter  S  defined  by 

r      _    kzdiL/bz  (20) 

v  -    u. 

Note   that    the   numerator  of  S   has    the    form 

kiLdu./az    3    If  3U./d2  (21) 

which,  according    to   the   Prandtl  model,  may   be    regarded   as 
representinn        tt^p,        U.^    as  here   defined    represents    that 
part   of      lt#     due   to    turbulent    fluctuations   associated   with 
the  mechanical   turbulent  mixing    length  jLfskz.     The    expression 
for   3  may    then  be  written 


v5 

The  drag  coefficient  defined  by  Cf  ^U-^/Vg    is  not  one 
that  can  be  measured  directly,  but  is  one  whose  value  can  be 
inferred  by  empirical  information  regarding  So   This  is 
dealt  with  in  more  detail  in  section  8„   This  interpretation 
of  U».e  as  kzdu./dz..  appears  to  have  been  given  first  by 
Busing er  [sj  . 

5.   Comparison  of  C/Ca  for  1953  and  1956  data. 

In  completing  Ca,  Lettau'3  empirical  formula,  eolation 
(17)  was  first  used  with  the  1956  data.   It  was  necessary 
to  compute  Z0,  using  ecuation  (9)  from  the  time-nearest 
neutral  wind  profiles  and  Vg  from  interpolated  surface 
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cherts.   In  computing  U.*and  C,  eouation  (8)  and  (12)  were 
used.   The  value  of  ftJ-140  ^s  also  Included  usinr-  equation 
(19).   The  results  for  all  data-days  employed  are  given  in 
table  1,  and  the  groupings  ol  ^-140  i-n  classes  are  given  in 
table  2.   The  graph  of  C/Ca  versus  Rii40  ls  given  in  fig.  3, 
with  Lettau's  curve  ^7J  superimposed.   The  agreement  with 
Lettau's  results  is  not  as  good  as  mi^ht  have  been  expected, 
For  instance  the  value  C/Ca-  1.75  was  found  it  Ri]_4Q-0 
Contrasted  with  Lettau's  value  C/Ca=lo00  at  Ki^oO'CL   More- 
over the  maximum  value  of  C/Ca  found  in  this  paper  exceeds 
Lettau's  maximum  by  a  comparable  factor  but  with  Ri]_40  - 
-0.055,  compared  with  Riioo  =  -0.007,  according  to  Lettau0 

Possible  sources  of  error  lie  in  the  values  of  Vg, 
which  affect  the  values  of  C  and  Ca  given  by  eouations  (12) 
and  (15)  respectively,  as  well  as  in  the  value  of  Z0»   The 
values  of  Z0  determined  here  ranged  from  0.35cm  to  5.1  cm, 
as  contrasted  with  values  quoted  by  Blackadar  et  al  [2]  of 
Z0~   0.6  cm  for  profiles  in  approximately  the  same  data<= 
period. 
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TABLE   I    -    Values   of  C/Ca  by    eauations    (3),    (12),    (15). 


Date   CST  Z0  (cms)  U*fcrn/sec)     Vf  (kts) 


10   Jul    1955 


C/Ca    Hl14o 


1305 

1.10 

24.67 

13.60 

0.036 

0.034 

1.08 

-Co45 

1405 

1.10 

35.85 

12.80 

0.055 

0.034 

1.65 

-0.17 

1505 

1.10 

31.24 

12.00 

0.052 

0.034 

1.52 

-0,28 

1530 

1.10 

»_-U  .  ^  J 

15.54 

0.040 

L  jOui. 

1.20 

-0.25 

1505 

1.10 

35.75 

11.  cc 

0.053 

0.034 

1.54 

=0.20 

1705 

1.10 

35.00 

10.80 

0.0^5 

0.035 

1.87 

-0.05 

1305 

1.10 

44.00 

10.52 

0.084 

0.035 

2.41 

-0  o  02 

2105 

1.10 

21.03 

12.00 

0.035 

0.034 

1.02 

0.06 

23   Jul    1956 


0930 

3.14 

49.26 

19.44 

0.051 

0.037 

1.38 

=  0.15 

1105 

3.14 

42.85 

19.00 

0.045 

0.037 

1.22 

=  0.21 

1205 

3.14 

46.69 

18.00 

0.052 

0.037 

1.39 

-0.16 

1305 

3.14 

52.92 

17.00 

0.0-^2 

0.038 

1.66 

-0.09 

1505 

3.14 

52.55 

16.20 

0.065 

0.038 

1.72 

=0.14 

1705 

3.14 

59.69 

16.00 

0.075 

0.038 

1.97 

-0o06 

1305 

3.14 

67.38 

15.80 

0.085 

0.038 

2.25 

=0.02 

2130 

3.14 

31.03 

9.52 

0.059 

0.040 

1.47 

0.04 

2305 

3.14 

31.86 

7.52 

0.085 

0.042 

2. 01 

0o03 

25   Jul    1955 


0205 

1.49 

2.15 

4.80 

0.009 

0.040 

0  o  22 

OolO 

0405 

1.49 

12.66 

9.60 

0.026 

0.037 

0.72 

0o05 

0805 

1.49 

50.49 

24.00 

0.042 

0.034 

1,22 

=  0.10 

("905 

1.49 

52.24 

26. OC 

0.040 

0.032 

1.24 

=0.01 

1230 

0.70 

58 .  37 

32.00 

0.038 

0.029 

1.27 

-0.02 

1405 

0.70 

65,24 

46.00 

0.028 

0.028 

1.02 

-0.02 

1605 

0.70 

52.80 

54.00 

0.023 

0.027 

0.35 

=  0.01 

2330 

0.70 

48.72 

44.80 

0.022 

0.028 

0.78 

Co  00 

6  Aup  1956 

1805 

0.35 

c,&  .  0«_, 

10.60 

0.056 

C.030 

1.84 

OoOC 

2105 

v»  vO 

11.65 

10.80 

0.022 

0.030 

0.71 

0.07 

2205 

0.35 

25.00 

11.00 

0.045 

0.030 

1.50 

0.01 
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TABLT'J  1    -    (Cont'd) 
Date   CST  Z0  (cms)    u£m/sec)     Vf  (kts)         C 

7    Auf    1955 


C/Ca    ftti40 


0405 

5.06 

57.72 

11.88 

0,097 

0.039 

2.47 

OoOl 

0805 

3.06 

63.14 

12.88 

0.098 

0.039 

2.52 

-0.04 

0905 

3.06 

109.10 

13.08 

0.167 

0.039 

4.50 

-0.02 

1105 

3.06 

69.66 

13.43 

0.103 

0.039 

2.67 

-0.08 

1530 

0.62 

59.79 

15.84 

0.075 

0.031 

2.44 

-0.04 

1605 

0.62 

64.83 

15.00 

0.086 

0.031 

2.77 

-0.01 

1805 

0.62 

36.02 

15.75 

0.046 

0.031 

1.47 

-0.01 

2205 

0.62 

23.87 

14.00 

0.034 

0.031 

1.09 

0.05 

2205 

0.62 

9.47 

12.40 

0.015 

0.032 

0.43 

0.13 

27    Auf    1956 


1230 

0.55 

42.96 

26.48 

0.032 

0,029 

1.12 

-0.07 

1505 

0.55 

38.09 

24.80 

C.031 

0.029 

1.06 

-0.02 

1705 

0.55 

40.08 

25.20 

0.032 

0.029 

1.10 

-0.02 

1930 

0.55 

1.79 

15.92 

0.002 

0,031 

0.07 

0.15 

2205 

0.55 

2.98 

26.60 

0.002 

0.029 

0.07 

0.04 

29  Aug  1956 

2135 

1.53 

24.09 

7.20 

0.067 

0.033 

1,75 

0.03 

TAELS    2     -   Class   means   of    Hi]_40  vs   C/C, 


Hi  140 

-0.43 

-0.37 

-0.31 

-0.25 

-0.19 

-0,12 

-0.06 

0.00 

0,06 

0,13 

c/ca 

1.08 

- 

1.52 

1,20 

1,45 

1.44 

2,04 

1.72 

1.11 

0.25 

12 


Pig.  5 

rtatio  of  geostrophic  drag  coefficient  for  non-adiabatic  ^na 
adiabatic  condi-tions  as  a  function  of  Richardson  number 
at  Z  = 140  cm,  at  O'Neill,  Nebraska  1956. 
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6.   Revised  formula  for  Ca 

Hitherto  Ca  has  been  computed  by  using  Lettau's  empirical 
function  for  Ca,  eouation  (15).   Since  Lettau  himself  realized 
that  the  numerical  constants  in  (15)  were  not  very  exact  |7j  , 
an  attempt  is  made  here  to  obtain  an  improved  relationship 
between  Ca  nnd  Ho  by  using  only  1956  O'Neill  data0   The  first 
3tep  again  consists  in  computing  Z0  for  the  time-nearest 
reutral  wind  profile  according  to  eouation  (9)0   As  before*, 
VX*o„  results  from  equation  (14)  with  k=0o38,  so  that 

U^s  0.061  au.*  (23) 

Eouation  (13)  is  then  applied  to  determine  Cao   The  results 

of  the  computations  are  shown  in  table  30 

An  empirical  function  for  Ca  was  then  obtained  by  the 

following  procedure: 

Let  lof  do  -  x  and  Ca^y  and  suppose  the  desired  function  has 

the  form 

y   =  be  (24) 

where  y  means  the  desired  function,  while  b  and  <g  are  cop- 
stants  to  be  determined.   Taking  the  lof  of  (24)  one  obtains 

Ky  =  1Ui>k»  (25) 

To  simplify  the  notation,  let  k)lsY  and  %*>bs.a,  ,    and  there- 
fore  the  desired  function  becomes 

Y  =  0-  +  CX.  (26) 

Using    the   least-squares  method,    the   normal    ecuations   are 

Y  -    OL  -CX.    ^    O 
?v  T  ?,  (27) 

Substituting  the  values  of  observations  (table  3)  in  the 
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normal  equations  one  obtains 

-3.51  -  a  -  7021c  =0 
-38L90-103.28a  -  785.18c -0 
from  which,  it  follows  that 

a  =  -0.05   or  b  =.  0,427 
Ca  -0.359 
Substituting  these  constants  in  (24),  and  further  by  chang- 
ing  y  to  Ca,  x  to  lof  Ro,    the  expression  for  Ca  c?n  be 
obtained 

Co.  =  0.427  e.*j>[-o.36g$og'Ro}        (23) 

The  fraph  of  equation  (28)  is  shown  in  fig,  4,  compared  with 
Lettau's  fraph.   Prom  fig.  4   it  is  evident  that  the  revised 
formula  for  Ca  tends  to  underestimate  the  draf  at  larger 
values  of  Rossby  number,  and  overestimate  it  at  smaller  Ro 
(lop  ito<7.1).   It  should  be  noted,  however,  that  there  was 
considerable  scatter  of  the  observations  relative  to  the  line 
of  best  fit. 

Jf  one  assumes  that  equation  (28)  is  a  universal  formula 
for   eostrophic  drag  coefficient  in  middle  and  higher  latitudes 
under  neutral  conditions,  one  can  then  write 

a^  -^JLi  =  %*  tvf-  (29) 

with  Ca  given  by  (23).   Hence,  adiabatic  wind  speed  at  level 


Z  is  obtainable  as  a  function  of  Vg,  Z0,  Z  and  latitude 
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Fig.  4 

Geostrophic  drag  coefficient  under  neutral  conditions  as 
a  function  of  surface  Rossby  number. 
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7.   A  revised  formula  for  C/Ca 

The  revised  formula  for  Ca  will  "be  used  in  this  section 
to  compute  C/ca,  still  using  the  previously  tabulated  values 
of  f eos trophic  winds  and  C  (table  3)„   The  results  are  shown 
in  table  4  and  the  resultin   values  cf  C/Ca  then  plotted 
versus  Ri]_4o  as  shown  in  fir  .  5.   The  procedure  of  fitting 
the  curve  is  shown  below. 

After  plotting  the  values  of  C/Ca,  a  further  restriction 
is  placed  on  the  empirical  function  C/Ca;  namely  that  C/Ca=l 
at  Hi=0.  I. oreover,  since  the  shapes  of  the  curves  for  C/Ca 
in  unstable  and  stable  conditions  are  different  (according 
to  the  scatter  cf  roints),  it  is  assumed  that  in  stable  condi- 
tions the  desired  function  has  the  form  j-z   eax  where  y^C/Ca 
and  x  a  ^i40»   ^n  unstable  condition,  y  is  assumed  to    cf 
the  form  y-   l4-k,X,6.  *    where  y=.C/Ca  and  x==Ri]_4os  a,  kj, 
k2#  are  constants  to  he  determined 0 

(i)  For  stable  conditions  (table  5),  with  Ri]_4o=x  and 
C/C  -  y,  if  the  assumed  function  has  the  form 

>=e  (so; 

then 

and  if  Ly*  Y 

Y  ■  =  0,* 

By    the  method    of   least-scuares,    the   normal    enuotion  will   be 

i*»  j-.i 

Solving    for  a, 

a  -  —  Z|.fS9 

IS 


Subs ti tut i rig    this    constant    in    (30)    and    chan ;inf   y   to   C/Ca 
x   to    Hi 140*    the    expression   for   C/Ca    Ls   obtained   as 


C/Ca=  exp[-^»ff9^3 


(32) 


19 


TABLE  b.  -  Computations  of  C/Ca  using  equation  (28)  for  Ca 


Date  CST 

log  Ro 

Ca 

c 

c/ca 

Rilli0 

10  July  1956 

1305 

6.799 

-2.51 

0.08I 

0.03b8 

0.0363 

i.oJ-5 

-0.U6 

lb05 

6.773 

-2.50 

0.082 

0.0350 

0.0560 

1.600 

-0.17 

i5o5 

6.715 

-2.U9 

0.083 

0.0360 

0.0521 

IMS 

-0.28 

1530 

6.860 

-2.53 

0.080 

0.03142 

0.0396 

1.160 

-0.25 

l6o* 

6.730 

-2.1,8 

O.08I1 

0.0359 

0.0530 

1.U75 

-0.20 

1705 

6.699 

-2.L7 

0.085 

O.0363 

0.06b8 

1.785 

-0.05 

1805 

6.687 

-2.1*6 

0.085 

O.0363 

0.0838 

2.320 

-0,02 

2105 

6.71*5 

-2,h9 

0.083 

0.0360 

0.0350 

0.975 

0.06 

23  July  1956 

0930 

6.500 

-2.  hO 

0.091 

0.0388 

0.0507 

1.305 

-0.15 

1105 

6.1)87 

-2.39 

0.092 

0.0393 

o.ob5i 

1.1U5 

-0.21 

1205 

6.h6b 

-2.38 

0.093 

0.0397 

0.0519 

1.305 

-0.16 

1305 

6.h39 

-2,37 

0.093 

0.0397 

0.0623 

1.570 

-0.09 

i6o5 

6.JU18 

-2.36 

0.09b 

O.ObOl 

0.o6b9 

1.615 

-0.1b 

1705 

6.1il3 

-2.36 

0.091 

O.ObOl 

0.07146 

I.860 

-0.06 

1805 

6.L07 

-2.36 

0.09b 

o.oboi 

0.0853 

2.150 

-0.02 

2130 

6.230 

-2.30 

0.100 

O.Ob27 

0.0591 

1.3^ 

o.ob 

2305 

6.08I 

-2.22 

0.109 

O.Ol|65 

0.0850 

I.830 

0.03 

25  July  1956 

0205 

6.215 

-2.29 

0.101 

0.0U31 

0.0090 

0.209 

OclO 

obo5 

6.516 

-2.1)0 

0.091 

0.0388 

0.026U 

O.680 

0.05 

o8o5 

6.721 

-2.1-7 

0.085 

0.0363 

0.0U21 

1.160 

-0.10 

0905 

6.9ii8 

-2.56 

0.077 

0.0329 

0.0b02 

1.220 

-0.01 

1230 

7.366 

-2.72 

0.066 

0.0282 

0.0368 

1.305 

-0.02 

ibo5 

7.523 

-2.77 

O.063 

0.0269 

0.028b 

1.055 

-0.0? 

l6o5 

7.593 

-2.80 

0.061 

0.0260 

0.0233 

0.895 

-0.01 

2330 

7.512 

-2.77 

O.063 

0.0269 

0.0218 

0.810 

OoOO 

6  Aug  19^6 

1805 

7.188 

-2.65 

0.071 

0.0303 

0.0557 

1*835 

0.00 

2105 

7.196 

-2.66 

0.070 

0.0299 

0.0216 

0.722 

0.07 

2205 

7.20li 

-2.66 

0.070 

0.0299 

o.ob5b 

1.515 

0.01 

7  Aug  1956 

Oho5 

6.295 

-2.32 

0.098 

O.OI4I8 

0.0970 

2.320 

0,01 

0805 

6.330 

-2.33 

0.097 

0.01415 

0.0979 

2.^60 

-O.Ob 

0905 

6.337 

-2.3U 

0.096 

O.OblO 

0.1666 

b«06o 

-0.02 

1105 

6.350 

-2.31 

0.096 

o.obio 

0.1032 

2.510 

-0,08 

1330 

7.111 

-2.62 

0.073 

0.0312 

0.0755 

2.b20 

-O.ob 

l6o5 

7.087 

-2.61 

0.07b 

0.0316 

0.086b 

2.730 

-0.01 

i8o5 

7.108 

-2.65 

0.071 

0.0306 

O.Ob37 

l.Sbo 

-0.01 

2005 

7.057 

-2.6c 

0.07b 

0.0316 

0.03U1 

1.069 

0.05 

2205 

7.000 

-2.56 

0.07b 

0.0329 

0.0153 

O.L65 

0.13 
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TABLE  lu   (Cont'd) 


Date  CST 

27  Aug  1956 
1230 

i5o5 

1705 
1930 
2205 

29  Aug  1956 
2135 


log  ho  .o.u9U^ 


7.390 
7.361 
7.368 
7.168 
7.392 


-2.73 
-2.72 
-2.72 
-2.6U 
-2.72 


0.065 
0.066 
0.066 
0.071 
0.066 


0.0278 
0.0282 
0.0282 
0.0303 
0.0282 


0.0325 
0.0307 
0.0318 
0.0023 
0.0022 


C/C, 


1.170 
1.090 
1.125 
0.076 
0.078 


hi 


1U0 


-0.07 
-0.02 

-=0.02 
0,15 

0.0U 


6.377       -2.35      0.095      0.0^06      O.O669      I.6U0        0.03 


TABLE  5.    -  C/Ca  in  stable  conditions 


L       Kr,4,*x.|        C/Ca-y6-       Jky;-Y? 


Sums 


M 


1 

0.057 

0.975 

-0.025 

-0.001 

0.0C3 

2 

0.036 

1.385 

0.326 

0.012 

0.001 

3. 

0.035 

1.830 

O.60I 

0.021 

0.001 

h. 

0.099 

0.209 

-1.565 

=0.155 

0.010 

5. 

0.051 

0.680 

-0.386 

-0.020 

0.003 

6. 

0.067 

0.722 

-0.326 

-0.022 

o.ooU 

7. 

O.Olii 

1.515 

0.I4I5 

0.006 

o.oco 

8. 

0.00P 

2.320 

0.812 

0.007 

0.000 

9. 

o.o5h 

1.069 

0.068 

0.00U 

0.003 

10. 

0.127 

0.165 

-0.766 

-0.01C 

0.016 

11. 

0.153 

0.076 

-0.271* 

-O.0li2 

0.023 

12. 

O.OiiO 

0.078 

-0.2U8 

-0.010 

0.002 

13. 

0.031 

1.6140 

0.1*95 

0.015 

0*001 

-0.283 


0.066 
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(ii)  For  unstable  conditions  (table  6),  the  scatter 
of  points  of  C/Ca  versus  Rii4Q  suggests  a  curve  of  the 
form 

y  =  i  4-  k.x.e^*  (33) 

where   -    fti]_4C  =  x  and    C/Ca=  y0      Again  C/Ca  has   been   forced 

to   approach  unity   as    Ri    approaches    zero0  If  we    let   y=lia*f 
and   then  take   natural   logs 

JL^  -  SL  k,  +  JL  x,  >  kax.  (34) 

jbf.   -Jfvwk.  ♦  k  x. 

Let    Am,  1  •=  Y      and     am,  k  s  a,  ,    and    substitute   in    (34) 

Y     =      A.   +    ka* 

The   normal   eauations   are 


Y-  a  -  ktx.  =  o 

WV  |W  M» 


(35) 


Solving  for  k]_  and  k2 

k1=  30.364  ,   k2  =  -13,324 
Substituting  these  constants  in  (23)  and  replacing  x  by 
•Ril40,   7  ty  C/Ca,  gives 

C/C^  =  *  -  so.a*/,*,*,  e>q>  [ts.sa/^J      (36) 

Note  that  C/Ca  of  equation  (36)  may  te  maximized  for  Rii^c 

r  |  ;  -0.075,  for  which  (0/Ca)n>ax   is  1.839.   This  is 
close  to  the  maximum  of  the  fitted  curve  of  fig.  5.   Lettau's 
graph  of  C/Ca  has  been  included  in  fig.  5  and  it  may  be  noted 
that  the  computations  of  this  3:aper  are  in  general  agreement 


with  those  of  Lettau,  although  his  data  did  not  span  as  lar^e 
a  ranfe  of  fti  values.   The  results  Indicste  that  C/noy   1  in  un- 
stable conditions,  whereas  C/Ca  <  1  in  stable  conditions0 
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TABLE  6  -  C/Ca  in  Unstable  Conditions 


I 

*•!=-**, 

y?=c/c^ 

!ft=tt-i 

fl/x; 

*f 

Y(x; 

1 

O.U63 

I.0I6 

o.ol*5 

0.097 

-2.333 

0.211* 

-1.080 

2 

0.165 

1.600 

0.600 

3.630 

1.289 

0.027 

Oc213 

3 

0.280 

1.1*1*5 

0.1*1*5 

1.590 

0.1*61* 

o,o7r 

0.130 

h 

0.251 

1.160 

0.160 

O.638 

-oJ>)*9 

0.063 

-0.113 

5 

0.199 

X.U75 

0.1*75 

2.390 

0.871 

o.ol*o 

0.173 

6 

0.053 

1.785 

0.785 

11*.  700 

2.688 

0.003 

0.11*2 

7 

0.017 

2.320 

1.320 

77.600 

a.  352 

0 

0.071* 

8 

0.ll*7 

1.305 

0.305 

2.080 

0.737 

0.022 

0.108 

9 

0.210 

1.11*5 

o.il*5 

0.690 

-0.371 

0.01*1* 

-0.078 

10 

0.15? 

1.305 

0.305 

1.9U0 

0.663 

0.025 

o.iol* 

11 

0.091 

1,570 

0.570 

62.600 

U.137 

0.008 

0.376 

12 

O.lljl 

1.615 

0.615 

a.  360 

1.1*72 

0.020 

0.208 

13 

o.o58 

1.860 

0.860 

114.800 

2.695 

0.003 

0.156 

111 

0.022 

2.150 

1.150 

<2.3x> 

3.957 

0 

0.087 

15 

0.103 

I.I60 

0.160 

1.550 

0.1*38 

0.011 

0.01*5 

16 

0.011 

1.220 

0.220 

20.000 

2.996 

0 

0.033 

17 

0.017 

1.305 

0.305 

17.950 

2.888 

0 

0.0l*9 

18 

0.019 

1.055 

0.055 

2.890 

1.061 

0 

0.020 

19 

0.0I42 

2.360 

1.360 

32.1*00 

3.1*78 

0.002 

Ooll*6 

20 

0.077 

2.510 

1.510 

19.600 

2.976 

0.006 

0.229 

21 

0.0L3 

2.1*20 

1.120 

33.000 

3.500 

0.002 

0.151 

22 

0.011 

2.730 

1.730 

157.300 

5.056 

0 

0.056 

23 

0.010 

l.5lo 

o.5lo 

5Ii.ooo 

3.989 

0 

O.OliO 

21* 

0.075 

1.170 

0.170 

22.700 

3.122 

0.006 

0.23U 

25 

0.021* 

1.090 

0.090 

3.750 

1.322 

0.001 

0.032 

26 

0.02k 

1.125 

0.125 

5.200 

1.6ltf 

0.001 

0.01*0 

iUXn      2.710  52.61*9      0.576        1.575 
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Pig.  5. 

Ratio  C/Ca  as  a  function  of  Richardson  number  at  Z = 140  cm, 
at  O'Neill,  Nebraska  1956  with  Ca  given  by  equation  (28). 


8.   The  normalized  — logarithmic  wind  sheai  S 

It  was  explained  in  section  4  that  if  one  knows 
S  sJ!l^i  ,  one  may  obtain  directly  Cf/C,  or  U^g/U^ 

LA-*-  ' 

Ellison  [2"]  has  fiven  a  prediction  equation  for  S  in  unstable 
cases 

S\.3gS*«i  <w) 

where  Z  is  height  arid  L  is  the  gradient  length  defined  as 

!»         u.«e  ao./di  (38) 

with  y'  defined  as 


Ki 


here  T  is  a  constant,  and  l<  and  i<m   are  eddy  diffusivity 


coefficients  for  heat  and  momentum  respeetively0   Divid 
both  sides  of  (37)  by  S4 


1  7 


-4 


-V'  a- 

Substituting  (20)  and  (18)  in  the  expression  for  •£-! 


no 


solving  for  S 

S  =  (-*'*•■)  (40) 

Elackadar  et  al  [_2J  suggest   the  value  TT  =  18  for  ui    ble 
conditions.   However  in  stable  conditions,  they  find  that 
equation  (40)  fives  no  better  approximation  than  t]  e 
Obukhcv  "lo£ +•  linear"  wind  profile  defined  by 

5  ~  \   +  ^!  Z  ,   <*  =consu^  (41) 

which  is  ordinarily  valid  for  Richardson  numbers  near  onec 
Apain,  it  may  be  verified  that  (41)  is  ecuivslent  to  an 
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ecuation 

5-  l  -hoC'fr  (42) 

Blackadar  et  si  [2]  have  tabulated  values  of  ex'1  at  O'Neill 

for  various  layers.   For  the  layer  1-2  m  under  stable  con- 
ditions, the  best-fitting  value  of  the  constant  was  oC  a 
which  will  be  used  for  the  stable  cases  in  this  paper0 

The  computations  for  S  are  shown  in  table  70   The  values 
of  5  are  presumed  applicable  at  the  height  assigned  to  Hi, 
that  is  at  Z  -  \Xrt\       or  Z-14104  cm0   The  curve  of  S  is 
shown  in  fif •  6.   The  result  indicates  that  S  < 1  in  unstable 
cases,  and  S >  1  in  stable  cases.   Whether  S  in  stable  condi- 
tions levels  off  with  a  further  increase  of  Hi,  is  not  cer- 
tain yet,  since  the  value  chosen  for  °0  may  not  be  valid  for 
larre  positive  values  of  Hi. 
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TABLE  7.  - 

-  Results 

of  computa 

tlons  fi 

Date  CST 

i 

K3-140 

Vg(kts) 

S(140 

10  July 

1956 

1505 

-0.46 

13.60 

0o572 

1405 

-0.17 

12.80 

0o708 

1505 

-0.28 

12.00 

0.638 

1550 

-0.25 

15.64 

0.652 

1605 

-0.20 

11.60 

0.683 

1705 

-0.05 

10.80 

0.850 

1805 

-0.02 

10.52 

0.935 

21C5 

0.06 

12.00 

1.515 

23  July 

1956 

0930 

-0.15 

19.44 

0.724 

1105 

-0.21 

19.00 

00676 

1205 

-0.16 

18.00 

0.715 

1305 

-0.09 

17o00 

0  784 

1605 

-0.14 

16.20 

0.729 

1705 

-0.06 

16.00 

Co836 

1805 

-0.02 

15.80 

0.916 

2130 

0.04 

9.52 

1.326 

2;305 

0.03 

7.52 

1.513 

25  July 

1956 

0205 

0.10 

4.80 

1.888 

0405 

0.05 

9.CC 

1.457 

0805 

-0.10 

24.00 

0.769 

0905 

-0.01 

26.00 

0.958 

1230 

-0.02 

32.00 

0.936 

1405 

-0.02 

46.00 

0.929 

1605 

-0.01 

54.00 

0.9<^0 

2330 

0.00 

44.60 

1.033 

6  Au£  19 

56 

1605 

0.00 

10.60 

0.994 

2105 

0.07 

10.80 

1.PC5 

2205 

0.01 

11.00 

1.123 

7  Aug  19 

56 

0405 

0.01 

11.88 

1.076 

0805 

-0,04 

12.88 

0.868 

0905 

-0.02 

13.08 

0.939 

1105 

-0.08 

13.48 

0.604 

1330 

-0.04 

15.  t  4 

0.867 

1605 

-0.01 

15.00 

0.958 

1605 

-0.01 

15.76 

0.960 

2005 

0.05 

14.00 

1.49C 

2205 

0.13 

12.40 

2.147 
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TABLE  7.  (Cont'd) 

Date  CST  Hi  140  v?  (kts)  s(140) 

27  Aur  1956 

1230  -0.07  26.48  0.808 

1505  -C.02  24.80  0.915 

1705  -0.02  25.20  0o914 

1930  0.15  15.92  2.375 

2205  0.04  26.60  1.362 

29  Aug  1956 

2135  0.03  7.20  1.278 
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Pig.    6. 

The   normalized    logarithmic  wind    shear  S  as   a   function   of 
Richardson  number  at   Z=  140   cm  based   on  1956   O'Neill, 
Nebraska   data. 
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9.   An  alternative  method  of  determination  of  U-* 

In  section  7,  C/Ca  was  computed  by  equations  (8),  (12), 

and  (28).   While  there  may  te  considerable  error  in  assuming 

a  functional  form  for  Ca  of  the  type  derived  in  section  6, 

there  is  also  a  possibility  of  error  in  the  computation  of 

U-*  by  equation  (8).   Hence  an  alternative  computation  of  U.* 

was  employed,  one  which  employs  the  expressions  for  S  of  the 

section  8.   This  approach  essentially  leads  to  a  value  of  U»* 

in  terms  of  the  wind-profile  Uts)  and  the  inter ral-mean  value 

of  S  in  a  sub-layer  of  the  surface  layer,,   Thus  for  example 

Q_  _^L  5*5:  -  Jl   dj*- 

and  integrating  from  Z  =  Z0  to  Z  -.  Z,  assuming  CL*  is  constant 
in  tbe  surface  layer  leads  to 


U  "  IT  5^1"  (42) 

where  S  is  the  intef  ral-mean  value  of  S  up  to  level  Zo 

As  in  the  previous  section  S  was  taken  in  accordance 
with 

S  =  \  *  ~"  l8&^*>)         unstable  cases 

(44) 
5-  I  +  9  &C*>  stable  cases 

Although  Priestly  [8]  has  raised  a  cuestion  regarding  the 
the  validity  of  a  constant  U.^  with  height  in  stable  condi~ 
tions,  it  was  assumed  that  a  surface  layer  at  the  depth  of 
at  least  2m  exists  for  both  stable  and  unstable  cases, 

Values  of  S  according  to  eouation  (44)  were  plotted 
using  Ri  (Z)  for  values  of  Z=  5,  35,  70,  140,  280cms0   The 
lowest  value  of   height  corresponds  to  a  level  at  the  ~eometric 


mean  of   the    rouphness   parameter   Z0   and    Z  =  25   cm,    usinp   the 

mown  values  of   2>u-/d2   and   £©/dz  for  this    layer0      The 
graphs   of   S   versus  %«*  2      was    plotted    and    the    interral-mean 
value   S    extracted.      Fir.    7    shows    an   example   of  an  unstable 
and    stable   case  which    occurred    at    1905   and    2105   CST,    respec- 
tively on   10    July   1956.      This   permits   one    to    solve   for  U^  i~ 
means    of 

U>=      \^°°\     .  (45) 

The  values  of  U-%  computed  "by  eouation  (45)  are  displayed  in 
table  8  alonf  with  the  values  of  <Xfes  listed  earlier  In 
table  1„  The  absolute  average  error  is  approximately  10^, 
although  the  alfebraic  error  is  Z%. 

Note  that  the  ratio  Cf/Ca  may  te  obtained  as  the  product. 

C^     C   Co. 

=  5  ^*/V3  (46) 

Co. 

The  value  of  S  as  a  function  of  Ri(Z)  may  be  obtained  from 
eouation  (44).   Values  of  U-*  may  be  obtained  from  eouation 
(45);  and  Ca  may  be  determined  from  eouation  (28) „ 

The  main  purpose  of  this  last  discussion  is  to  empha- 
size that  Cf  is  not  eoual  to  Ca.   Listed  below  in  table  8 
are  some  computed  values  of  Cf/Ca<, 
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Fig.  7 

The  normalized  logarithmic  wind  shear  S  as  a  function  of 
£*-z  for  unstable  and  stable  conditions  for  10  July  1956, 


TABLE  8.    Cf/Ca  in  unstable  and  stable  cases 


Date   CST 


^140 


U£:; /sec)      (Cf/Ca)(  6i^/sec)froni ( 3 ) 


10  July 

1956 

1505 

-0.46 

0.572 

2501 

0.605 

24  0  7 

1405 

-0.17 

0.90C 

35.0 

1.106 

1 R  C 

.505 

-0.28 

0.865 

«^X  a  w 

0.925 

31 0  2 

1605 

-C.20 

C.9C0 

51.6 

1.038 

35.8 

1705 

-0.05 

0.950 

34  „  8 

1.515 

35  oC 

1805 

-0.02 

0.970 

43.0 

20100 

44o0 

23  05 

0.06 

1.150 

20o5 

1.435 

21.0 

25  July 

1956 

1105 

-0.21 

0o850 

«^D  .  ' 

00648 

42.9 

1205 

-C.16 

0.880 

38  0  7 

0o773 

46  0  7 

1505 

-0.09 

0.885 

4403 

1.025 

5209 

1605 

-0.24 

0.885 

42.6 

0.956 

5205 

1705 

-0.06 

0.950 

46.7 

1.215 

59.7 

2505 

0,04 

1.115 

27.1 

2.030 

<uL  0  8 

25  July 

1955 

0405 

0.05 

1.100 

1609 

1.325 

12.7 

5  Aug  1956 

2105        0.07    1.200     12.1      1.210 

7  Auf  1955 

2005        0.05    1.100     22.6     1.517 


11.7 


23    ° 


27    Aug    1956 
2205 


0.04 


1.00 


10.6 


0.386 


7.93 
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10.   Summary  and  conclusions., 

A  functional  relationship  between  C/Ca  and  Ri  was 
obtained  by  usinf  an  expression  for  Ca  suf rested  by  Lettau0 
Better  agreement  was  obtained  after  makinr  some  revisions 
in  Lettau's  formula  for  Ca<>   An  expression  for  Ca  has  the 
practical  significance  that  one  can  yet  directly  the  ratio 
of  a  windspeed  at  some  anemometer  height  to  a  surface  peo- 
strophic  wind  by  knowinr  the  rou.'hness  parameterQ 

By  usinf  the  revised  expression  for  Ca,  a  new  functional 
relationship  between  C/Ca  and  Hi  was  obtainedo   This  function 
has  different  analytical  forms  in  unstable  and  stable  condi- 
tions. 

Another  way  of  obtaining  C/Ca  was  based  upon  computa- 
tion of  U^.  usinp  Blackadar's  results0    By  this  means  a 
wind-profile  expression  was  obtained  which  varies  with 
stability.   Prom  this  new  procedure   Dne  can  :~et  the  value 
of  Cf/Ca  without  using  C/Ca. 
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APPE1IDIX 
Computation    for  Ca   using    Lettau's    formula 

The   original    ecuation    is 


C<x     ^ 


By    cross  multiplying 

or  C  2.-24  +0  5o4/cA) 

O   iOii/Cc 

=  73.5   A   10 

then 

S"-*0  -  exp  ( f J 

-   exjo  (  G^tpVCo.) 

Let   0.2395/Ca  =  y   and    substituting    in   the  atove    ec-ation 

o, 3330-/173. gRo  -  V&V.,  o 

Since    y  *  O 

Taking  the  lors  to  the  base  10 


(47) 


J^(5^«.)=  Lj  J  +  0.^.3/,3j/  (48) 

Let    lor  f^-SrRo)—      l<     ■=  constant,  and    replace    lo(    y  by 

0o4?42J^y      and    substitute    in    (48) 

f(y)=  o^Ail^y  +  0.43^3^-  K-©  (49) 

To  pet  the  root  of  ecuation  (49),  we  will  use  the 
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Newton- rtephson  method 

M**y.--  -ft—  (so 

Since  t  <V>  =  0.  4343/y -»-  00  4343,  therefore  by  substituting  this 
in  (50)  we  obtain 

Substituting  lop  (°'yj3|*  *°)         for  K  and  divide  the  fraction 

by  0.4345  , 
0.4343 

»/y,   +  i  

Chancing    **\y4-    to    2.3026    lop   y^_ 

*,3c*6  L3y(  +y,  +  6.S317  -  *,3oxs-lg  i  c 

yU*    -    35'    -  l/y.    +|  (51; 


37 


BIBLIOGRAPHY 


1.  Barad,  M.L,  Geophysical  research  papers  no.  59 j  vols  1 
and  2,   GRD  1958. 

2.  Blackadar,  A.  K.,  et  al,  Structure  of  turbulence  and  mean 
wind  profiles  within  the  atmospheric  boundary  layer, 
Pennsylvania  State  University,  1960. 

3.  Busin^er,  J.  A.,  On  the  structure  of  atmospheric  surface 
layer,  J.  of  Meteor.,  vol  ]2,  number  6,  553-561,  1955=, 

4.  Davidson,  B. ,  and  M«  L.  Barad,  Some  comments  on  the 
Deacon  wind  profile,  Trans.  Am.  Geophys.  Union,  19560 

5.  Haltiner,  G.  J.,  and  F.  L.  Martin,  Dynamical  and  physical 
meteorology,  McGraw-Hill  Book  Company,  inc.  1957 . 

6.  Lettau,  H.  and  B.  Davidson,  Exploring  the  atmosphere's 
first  mile,  vol.1,  Pergamon  Press,  New  York,  1957. 

7.  Lettau,  H.,  Advances  in  geophysics  no.  6,  241-257, 
Academic  Press,  New  York,  1959. 

8.  Priestly,  C.  H.  B.,  A  determinant  hypothesis  for  the  super- 
adiabatic  wind  and  temperature  profiles,  Quarterly  Journal 
of  the  Royal  Meteorological  Society,  vol.  86,  No.  368, 
232-236,  April  1960. 


33 


